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An algorithm is described for sparsely updating the filter coeffi- 
cients of an adaptive filter of an adaptive audio system to facilitate a 
high sampling rate without excessive processing capacity. Computer si- 
mulations are used to illustrate the potential for using an inverse filter 
matrix to manipulate the position of stereophonic images. The possibil- 
ities suggested could be used in conjunction with existing stereophonic 
recordings in order to improve stereophonic imaging in situations 
where the relative locations of loudspeakers* and listener are con- 
strained to be non-ideal. In one arrangement, real sound sources (SI, 
S2) provide to the listener the impression that the sources are at more 
widely-spaced positions (VI, V2). In another arrangement, employing 
four sound sources (SI, S2, S3, S4) the impression is created at two lis- 
tener locations (Ml, M2; M3, M4) that the sound is coming from a 
pair of sound sources positioned syomietrically of those listener loca- 
tions. 
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ADAPTIVE AUDIO SYSTEMS AND SOUND REPRODUCTION SYSTEMS 

This invention relates to adaptive audio systems, and to sound 
reproduction systems incorporating multi-channel signal processing 
techniques. 

SUMMARY OF THE INVENTION 

Broadly speaking, the first three aspects of the invention are 
concerned with sound reproduction systems arranged to generate 
'virtual source locations' of sound at positions other than those of 
actual loudspeakers enployed to reproduce a sound field. 

According to one aspect of the present invention a sound 
reproduction system comprises a plurality of loudspeakers which are 
arranged asymmetrically with respect to a listener location, the 
loudspeakers being driven through a filter means by a plurality of 
channels of a sound recording, the filter characteristics of the 
filter means being so chosen as to create at the listener location a 
local sound field which is substantially equivalent tp the local 
field that would result from playing the sound recording through a 
plurality of loudspeakers driven without filters and positioned at 
virtual source locations that are substantially symmetrically 
positioned with respect to the listener location. 

Thus, the filter means is arranged to compensate for the asymmetric 
positioning of the loudspeakers with respect to the listener 
location. 

For example, with reference to Figure 1, a transfer function matrix 
C(2) of electro-acoustic transfer functions relates loudspeaker 
inputs to the outputs of microphones placed at the location of the 
listener ears in the sound field. A matrix H(2) of inverse digital 
filters can be used to process conventional two-channel stereophonic 
recorded signals prior to transmission by the loudspeakers in order 
to produce the desired effect, as specified by a matrix A(z) at 
microphones placed in a sound field. Matrix A(z) is the transfer 
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function matrix relating the recorded signals to the signals that 
are desired to be reproduced at the microphones, and the matrix A(z) 
in accordance with the first aspect of the invention is selected by 
assuming that we have the said specified virtual source locations. 

According to a second aspect of the present invention a sound 
reproduction system comprises at least two loudspeakers which are 
driven through a filter means by at least two channels of a sound 
recording, the filter characteristics of the filter means being so 
chosen as to create at the listener location a local sound field 
which is substantially equivalent to the local field that would 
result from playing the sound recording through the loudspeakers 
driven without filters and positioned at virtual source locations 
that are more widely spaced from one another than the actual spacing 
of the loudspeakers. 

Thus the filter means is arranged to create an impression at the 
listener location that the loudspeakers are more widely spaced than 
is actually the case. 

According to a third aspect of the present invention a sound 
reproduction system comprises four loudspeakers which are arranged 
at spaced-apart positions to create a sound field, first and second 
predetermined listener locations within the sound field, the 
loudspeakers being driven through filter means by two channels of a 
sound recording, the filter chararteristlcs being so chosen as to 
create at both the first and second listener locations a respective 
local portion of the sound field which is substantially the same as 
the sound field portion that would be produced at that respective 
location by playing the unflltered channels of the sound recording 
through a pair of loudspeakers positioned synmetrically with respect 
to the respective location. 

A fourth aspect of the invention is concerned with an adaptive audio 
system and to a method of updating the filter coefficients of the 
adaptive filter of the system. 
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According to the fourth aspect of the present invention an adaptive 
audio system comprises an adaptive filter having a plurality of 
alterable filter coefficients, a processor implementing an algorithm 
which cyclically performs in turn a filtering operation utilising 
the-filter, and a filter updating operation to update the filter 
coefficients, in which the algorithm is arranged in a cycle thereof 
to adjust (if necessary) only a limited number of the filter 
coefficients before performing a filtering operation. 

This aspect of the invention enables the use of a filter with a 
relatively large number of coefficients whilst facilitating a high 
sampling rate. 

Preferably only one coefficient of the filter is adjusted in each 
cycle. 

The algorithm is preferably the LMS algorithm [2] or the filtered-x 
LMS algorithm [3]. In the latter case, the filtered reference 
signal employed can be calculated on a single-tap basis. 

It is shown hereafter that the sparse update implementation method 
of the fourth aspect of the invention reduces the operation count by 
a factor of 2 in the case of the LMS algorithm and by a factor of 3 
in the case of the filtered-x LMS algorithm. More importantly, the 
bulk of the computational work that remains to be done is mainly 
related to the actual filtering operation. These calculations can 
be performed by a dedicated filtering unit external to the main 
processor and then most of the processing time can be dedicated to 
update the filter. This technique can be easily extended to complex 
systems based on a number of adaptive filters having arbitrary 
lengths. This is an ability which is highly desirable in 
multl -channel applications, such as those of the first, second or 
third aspects of the invention. 
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INTRODUCTION 

the sound PJ^'y;^ " * P response- in time domain terms, an impulse applied 
response and » (at a little time later) at the 

this pnnaP^e ^ the use of multi-channel inverse filters in the 

considerable potential ° badceround to the Hlter design problem 

reproduction of sound, in multi-channel problem is 
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sound field in order to compare the reproduced signals w, J the signals that aio 
Sd To be reproduced. The fact that the filte« can a«refore be d«^^^^ 
aSptively opens up the possibility of tailonng the individual filters to Ae 
req^remLts of a particular listening space However, the filters do not 
nSSsarilv have to be adaptive and modified to accommodate changes in 
~om a^mtics; the design process can equajly well be mord^ to 
^^ish fixed filters uied, for example to modify the position of stereophonic 
fm?^?ta in<« SiteSninent sys^. It U the latter topic that we address m 
SSfnT wheiTSom^r simulaiions are presented which demonstrate some 
pos'biS. in particular, it is shown that * -^^.i'j--^^^;, 
usid to provide "virtual source" locations at positions other than those of 
Z J^noudspS3ce«. Furthermore, the possibUity is examined of providing 
xJ^lS? •ide^^SeJine positions"; it is demonstrated that an appropriately 
multiple Ideal P?*^ operate on the two channels of a 

Z^Jin^e reproduced field. The ultimate test of these possibiUties wiU o^ 
Sirfe Se psySoacoustical. In the meantime however, attenuon will be 
concentrated on the physical aspects of the problem. 

1. MULTI-CHANNEL INVERSE FILTERING USING A LEAST SQUARES 
FORMULATION 

LI Background to least squares filter design 

The desiitn of digital filters for single channel equalisation is most readily 
aooroaS usini traditional "least squares" methods. This technique has its 
Ste Si the Ssfcd approach of Wiener HI in which the impulse response of a 
mer S *nstS Tbe causal and designed in -f-.J J^-^/^^^^^^^ 
averaged squared error between the filter output and the desired ^^jP"^ 
Xemine equations which have to be satisfied to ensure an optim^ design 

I^'SrT?aSle when working In ^^^^^ ^^^'^^ .^^^Smte^l^ 
^ least sauares method has become a standard techmque m digital filter design^ 
t*^o'" *e lliS algorithm of Widrow Hoff Ul provid^ an^^^^^^ 
numerical technique for rapidly adapting the coeffiaents of 
orovide the optimal impulse response. In acoustics, this algorithm "q""*? * 
ft^CT modification before it caJ be utilised. In the loudspeaker equahsaUon 

indeoendentlv for use in feedforward control by Widrow (5) and 
«;!^or?oiid by Burgess 161. In many problems invo vmg the active conwl 
of sound and vibration it is often necessary to use ""^^^e ^pu" 'l^"^^, 
that the control filters are designed to ensure the 

appropriate measure of error at multiple points m space (71. J^qmrem^ 
l£l EUiott and Nelson 181 to generalise the filtered-x "^/8°"trthe MuS 
multiple errors. The resulting algorithm has become ^.^ '•'^^^'llv^^^^ 

SJor LMS al^oritiim 191 and it has been extensively utdised m a variety of 
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appUcations involving the active control of sound and vibration (see Nelson and 
Elliott 110] for a full accounO- 

The Multiple Error LMS algorithm has been stiU further generalised 
sD«ificallv appUed to problems in audio system equalisation by Nelson a al fll, 
^zTlirAat SSrk t^^ formulation of the problem was generalised to incorporate 
myplelnUrsignaL. This is the case found in stereophonic reproduction for 
e^vi wCe Ae algorithm has since been applied with considerable succe» 
TTtJ; gei^al formulation of the least squares filter design P^^tem ^11 Jgam 
be presented here before describing further potential applications of the 
technique to the reproduction of sound. 

■L2 Solution of the multi-channel problem 

The multi-channel problem is illustrated in block diagram form in Rguje J 
WorlTg irScrete time, we have K recorded signals .r,(n) comprising the 
ve^tTxV) These are transmitted via M loudspeaker channels whose inputs are 
«^en by v„(n) whii are the elements of the vector y(n). The resulting signals 
fre tran'sStVd via a matrix C(z, of ^^-troacoustic trans er ^^cti^^^ 
detected by L microphones whose outputs are given by z^M. these bemg tne 
elements of the vector zM. We introduce a K x M matrix H(z) of FIR digital 
S which oLates on the K recorded signals prior to 
foX"^e« T^ie coefficients of the fUters in the rnatrix are d«'g"f " ^ 
mSse the weighted sum of the mean squared values of the error signals 
S„) -^e ^tS signal is defined as the difference between the desired signal 
ff„) a?d the reproduLi signal The desired signals rf,(n <<:ompnsing die 
vector din)) are in turn specified by passing the recorded signals x,in) through a 
KxL matrix A(-) of filters. The filter matrix A(2) thus specifies the desired 
fi^L wSlsl l^Lt the most general method of ^Pl^^t^^^^^^^^ 
thfelements of A(r) wiU in general include W'j^P"** ^^S^^ ^^LSed 
that the desired signals are in some sense delayed vers««« of ti^ «oomea 
signals. This is dilrly necessary if reductions in mean «l""!i«"?r ^ 
achieved when the elements of CM are ^<>^-«'';;;"';;";P^'^^ ''^^ 
which A(z) may be specified will become clearer m the next section. 

zj(n) can be expressed as 



Where h„M is the i'th coefficient of the HR filter whose input is the k th 
recorded s^gjal and whose output is the input to the m'th lou6s^^ ^ 
Ster is asstSned to have an impulse response of / samples m duration. In vector 
notation we can write equation (1) as 
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(5) 



(6) 



- h1-(o) r,(n) * hT(l) r,(« - 1) - - • hTtf - 1) t, (n - / + 1) (2) 

where we have defined a composite tap weight vector and a reference signal 
vector respectively by 

hr(0 = [fcn(0 "21(0 - . /.Mi(0 Ihii (0 /.22(0 • • h^iO \h,K^O H,x(0 . . hMic(0] (3) 

r,T(n) = [rin(«) ''m («> • • nm\M I r|2(n) ''/22<n> • • 

• • nM2 (") I n\K ("> ''I2K • • nMK(«)i 
A further composite tap weight vector can be defined which consists of all the / 
tap weights of aU the L X M filters. Thus is given by 

wT=[h'"(o)h'"a)..hT(/-l)] 
The L'th order vector of error signals can now be written as 

e(n) = R(n)w - d(M) 

where the matrix R(n) of filtered reference signals is given by 

nT(n)r,T(«-l).. riT(H-f + l) 
r2T(n) rj^n - 1) . • • rj^Cn - / + D 

Tj(n)Tinn-\)...Tjin-l*l) . 

With the vector of error sisals defined ^jf J^f ;,Tghrvi^^^^ 
determine the optimal value of the com^ite «PJ«^^*Ju^~ the 

S« c««ld«.B. Thus w. dtf« . COS. fa«aion g.»«. by 
I . Sltnn) W, e(«) ♦ ••'•(«) W, wW] 
wh« W. ^ W. .r. (gen«ay diagonal) wtighdng J,^^^ 

Sei^isr cn^r^T^rJ-w-'^r^i.r's.ui Luring w, . 

symmetric, / can be written as 

/ = E[dT(„) W, d(n) - 2 wTRT(„) W, d(n> * wWn) W,R(«) * W„)wl (9) 



R(n) = 
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The gradient of / with respect to w, also assuming that W^, is symmetric, can be 
written as 

^=:2E[(Rnn) WeRCn) + W - R^(n) W, d(n)] (10) 

Thus the solution that ensures that 3//3w is zero is given by the optimal tap 
weight vector 

w„ = [j^RHn) We R(n) + W^]}"' {E[RT(n) W, d(n)l} (H) 
The corresponding minimum value of / is given by 

Jo = E[dTM We d(n)] - E[dT(«) W, R(n)]w, " ^^^^ 

The optimal tap weight vector can clearly be found by inversion of the matrix 
EfRl'dt) WeR(n) + W„] which must be positive definite for a unique minimum 
to exUt. In the case when W^ = I (the identity matrix) and W„ = 0, this matrix 
has a block Toeplitz structure and efficient numerical schemes exist for »ts 
inversion (14). The other approach is to use the Multiple Error LMS algorithm. 
This has its origin in the method of steepest descent in which the minimum ot 
the function is found iteratively by updating the coefficient vector w by an 
amount proportional to the negative of the gradient of the function. First note 
that using equation (6) in equation (10) allows the gradient to be written as 

^ » 2E[RT(n) We e(n) + W^w] 

Following Widrow and Hoff [21, we now make the assumption that the filter 
coefficients are updated by an instantaneous estimate of this gradient, which is 
^en by dropping the expectation operator in equation (13). Thus the tap weight 
update equation becomes 

w(n + 1) = w(n) - a W„w(n) - o K^n) Vie «(«> 

where a is a convergence coefficient. Equation (14) thus spedfies a stople ajd 
readily implemented algorithm for iteratively converging solution for toe 
optimL aSffident vecS,r. As pointed out by Elliott et a 
-effort- weighting W„ can be shown to be equivalent to leaky LMS 
algorithm in which case in the absence of an error term e n), tiie coeffiawit 
vSor w would decay away to zero. Note that the implementation of he 
alRoritiim requires the generation of the filtered reference signals r,^(n) which 
comprise the elements of the matrix R(n). These can be gen^ted by passing *e 
recorded signals Zk(it) through HR filters which give an estimate of the transfer 
function C,!,(z). These in turn can be identified by using a broadband training 
signal passed through the m th loudspeaker to the / th microphone with the 
LMS algorithm used to adapt the filter coefficients. 
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13 Relationship of the least squares solution to methods for finding exact 
inverse filtexs 

An alternative approach to inverse filtering in room acoustics is that proposed by 
Myoshi and Kaneda (151 in the forn\ of the Multiple-Input/Output Inverse 
Rlteririg Theorem (MINT). In that work, it is demonstrated that a pair of filters 
can b4 designed which can be used on the inputs to two loudspeakers which are 
both used to traiOTiit a given recorded signal to a specific point in space. The 
filters can be designed to ensure that "perfect" equalisation of the transmission 
path is produced, even when the transmission paths between the loudspeaker 
inputs and the point where equalisation is required are non-minimum phase 
(after an appropriate bulk delay has been subtracted; a point which is not 
altogether clear from Myoshi and Kaneda s paper [15J). An analysis of the 
relationship between MINT and the least squares approach presented above has 
been presented by Nelson et al (16). This leads to a useful result which may have 
some bearing on the choice of the values of fC, M and L in a given application, 
together with number of coefficients /.used in the inverse filters. 

It is assumed at the outset that the transmission paths C/,„(2) can be adequately 
represented by HR filters having / coefficients. It can then be shown that to 
produce exact equalisation of the transmission channels trom K recorded signals 
to L = K points in space requires that the number of coefficients in the inverse 
filters is given by 



(15) 



The full derivation leading to this result is presented in reference (16). However, 
it should again be emphasised that the choice of 1 given by equation (15) ensures 
the exact equalisation of the J-coeffident FIR filters representing the transmission 
paths C/^(2); equalisation of the real transmission paths therefore assumes that 
ail these paths are exactly represented by /-coefficient filters. Nevertheless, the 
analysis presented in reference {16J gives some indication of the number of 
coefficients / required in the inverse filters and demonstrates furthermore that in 
order to realise an exact inverse, it is required that M > L (i.e. a greater number of 
loudspeakers is required than the number of points at which equalisation is 
attempted). ThU result is therefore consistent with the work of Myoshi and 
Kaneda (15) who show, for example, in the case M = 2; L = 1, that / = (/ - 1). The 
work presented in reference (16) generalises Myoshi and Kaneda's result and also 
suggests that the Multiple Error LMS algorithm can be used to find the required 
solution for the coefficients of the inverse filters. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram for the analysis of multi-channel 
inverse filtering problems in sound reproduction, 

Figure 2 shows a geometrical arrangement of real sources (SI, S2) 
whose outputs are prefiltered in accordance with the invention in 
order to produce signals at microphones (Ml, M2) which appear to 
originate from virtual sources (VI, V2)- 

Figure 3 shows the impulse responses from inputs (dashed line) 
and (solid line) to the two microphones shown in Figure 2 before 
(Control = OFF) and after (Control = ON) the introduction of the 
filter matrix H used to pre-filter the outputs of sources SI and S2. 
The relative arrival times of the impulses at the microphones are as 
if they originated from virtual sources VI and V2. A modelling 
delay of 6.25 ms included during the simulations has been removed 
from the prefiltered responses (lower plots). 

Figure 4 shows the impulse responses of the elements of the filter 
matrix H designed to pre-filter the outputs of sources SI and 52 
shown in Figure 2. The dashed lines show the impulse responses of 
the filters operating on input x^ and the solid lines show the 
impulse responses of the filters operating on input Xg. 

Figure 5 shows frequency response functions associated with the 
composite system H(z)C(z) whose impulse response is illustrated in 
the lower traces of Figure 3. The phase responses have been 
calculated relative to that associated with the first arrival in the 
time domain. Note the small difference in magnitude response and 
the substantial difference in phase response corresponds to that 
which would be expected to be produced by the virtual source 
locations shown in Figure 2. 
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Figure 6 shows a further geometrical arrangement of real sources 
(SI, S2) whose outputs are preflltered in accordance with the 
invention in order to produce signals at microphones (Ml, M2) which 
appear to originate from virtual sources (VI, V2). 

Figure 7 shows impulse responses from inputs (dashed line) and 
(solid line) to the two microphones shown in Figure 6 before 
(Control = OFF) and after (Control = ON) the introduction of the 
filter matrix H used to pre-filter the outputs of sources SI and S2. 
Note the increased time between the arrivals of the impulses 
produced by the location of the virtual sources VI and V2. 

Figure 8 shows impulse responses of the elements of the filter 
matrix H designed to pre-filter the outputs of sources SI and S2 
shown in Figure 6. The dashed lines show the impulse responses of 
the filters operating on input x^ and the solid lines show the 
impulse response of the filters operating on input Xg. 

Figure 9 shows frequency response functions associated with the 
composite system H(2)C(2) whose impulse response is illustrated in 
the lower traces of Figure 7. The phase responses have been 
calculated relative to that associated with the first arrival in the 
time domain. Note the small difference in magnitude response and 
the substantial difference in phase response corresponds to that 
which would be expected to be produced by the virtual source 
locations shown in Figure 6. 

Figure 10 shows a source and microphone layout used in simulating 
the production of multiple stereo images. The two signals x^ and 
are filtered in accordance with the invention with a matrix H{2) 
prior to transmission via sources SI - S4. 

Figure 11 shows impulse responses of the composite system H(2)C(z) 
for inputs x^ (solid line) and X2 (dashed line) when H(2) is 
designed to ensure that x^ is reproduced at microphones Ml and M3 
and that is reproduced at microphones M2 and M4 (see Figure 10). 
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Flqure 12 shows frequency response functions corresponding to the 
impulse responses of Figure 11. 

Figure 13 shows the impulse response of the composite system 
H(z)C(z) when (solid line) and (dashed line) are filtered by a 
2x2 matrix H(z) prior to transmission via sources SI and S2 in 
Figure 10. H(2) is designed to ensure that x^ and X2 are reproduced 
at Ml and M2 respectively. Note the degradation of the response at 
M3 and M4. 

Figure 14 shows frequency response functions corresponding to the 
impulse responses of Figure 13. 

Figure 15 shows impulse responses of the composite system H(z)C(z) 
when x^ (solid line) and Xg (dashed line) are filtered by a 4 x 2 
matrix H(z) prior to transmission via all four sources in Figure 10. 
H(z) is designed to ensure that x^ and are reproduced at Ml and 
M2 respectively. 

Figure 16 shows frequency response functions corresponding to the 
Impulse responses of Figure 15. Note the improvement in cross talk 
cancellation at Ml and M2 compared to the results of Figure 14. 

Figure 17 shows at (a) a block diagram for the adaptive system 
identification problem. The filter G must be adapted so that its 
output approximates the output from the plant, and at (b) a block 
diagram for the adaptive system equalization problem. The filter H 
must be adapted so that the output from the plant approximates the 
output of a given target system which, in general, includes a 
modelling delay that ensures the existence of a causal optimal 
filter. 

Figure 18 shows the frequency response function of a loudspeaker in 
an anechoic chamber. Unequalized (dashed line) and equalized using 
in accordance with the invention the sparse update adaptive 
algorithms (solid line). 
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2. POTENTIAL APPLICATIONS OF MULTI-CHANNEL TECHNIQUES IN 
THE IMPROVEMENT OF SOUND REPRODUCTION SYSTEMS 

2.1 The use of multiple eirois in single channel response equalisation 

Of course the simplest equalUation problem corresponds '° ^^^J^^^Jl^V^, 
inverse filter is iitroduced on the input to a 1°"^^?^ '^"^^^ 
designed to ensure that the signal reproduced at a point in SP*« ^^^JJ^f " 
Sle (in the least squares stnse) to a delayed version of *e mput signal. A 
c5r««i* *.^amnle of such an approach is that presented by Kunyama and 
U^J who ie Ae sin^^channei filtered^ LMS algorithm to design a 
S^rtte iSation of the on-axis response of a three-driver louckpeaker 
SJ^andIi^crossK>ver networks). A 512 coefHdent HR filter operating at a 
rTte oTsf kHz succeeded in pRKiudng a flat «nplihKie and linear pha« 
^ase £S.een 200 Hz and 12 kHz. No results r'^^^.^P^'^^Ji^S^^J^ 
^Lements of the response of the system at spatial positioiis other than those 
JT^nSiSsation was achieved. This is a cruda^^ issue m response 
eauaEo? problems; the effect of the equalising Hlter on the global 
of the system. This problem has been addressed ^ Wikon 1181 m 
KoSSrofloudsi^ equalisation. It was demonstrated that imnimi«jig a 
w^^^tmi of ^red erSrs derived from measurements of both on- and off- 
SSSonses was successhil in improving the off-axis response also, albeit at the 
SSSTSe^p^lements in oVi^ response. Similar comments apply to 
S^Jspo^^SaSsation. Famsworth «/ (19), for e«unp^e used compuj^ 
sSSa^ of an endosed sound field to investigate the effect f ^J^J^^^g *^ 
Tnsmission response to one point in a room on ^P^J^** «*5 
the room. The results predicted a potentially severe degradation of 'esP"«« 
outsJ^ zoai surrounding the equalisation points. The subjective respor«e to 
Ss Sect however, has «t to be quantified. An attempt to ensure global 
Suaita?on" w^^^ by Elliott «d Nelson 1201 who also "^eci comp^^r 

S^ulations to study the equalisation of the low frequency acjushc heW in a 
STerior. In that cie, muldple errors were again 

«;«*«^ie Koincr cnprified at a number of points m the sound held. These aesirea 
S tS?e^SS^S^'e«to^f the input signal xvith the delays di««n to be 
:?^to"S^ iSc propagation delays^om the ^^^'^^^^^^^ 
to the points at which the desired signals were speafied. This ^^^^^l^'^^ 
^l^alequalisation- in the sense that the frequency resjHJnse was »mpro«^K^ at 
Jl the 4SSsSon points, but of course was not made perfect at any one point. 
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23 Response «il«alis»llo»'»««'«>P''»™"''"°''°'*'°" 

tofter th«, ««n«p«»g (the umeatoble) go^ of Perfect f'^""' 'fl'"'^,^;; °2 

Sly the recorded signal from the left *annel and »he nght ear ^^^^ F^r« 
vice-versa, in addition to equalling ';*i^"*"f ^ "^^^ 
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capabilities of the original analogue implementation beikg graphically described 
bv Schroeder (22). Sudi a system is wdl suited to the reproduction of bmaurally 
recorded sound fields, such as those recorded with the use of an artificial head. 

However, most modem recordings do not use such technique and tfie typical 
twMAankel stereophonic recording is a carefuUy mixed amalgam of mdividual 
signals, with left and right channels attributed with signal components m order 
to ma^mise the effective stereophonic Ulusion, largely m a<«>';d«« wtti Ae 
subjective judgements of the producer of the recording. If one accepts that sudi 
rSng tech^ues wUl not readUy be dianged, the potential for the further 
2SSon of equalisation tedmiques becomes focussed on providmg re^^ 
l^rover^iS in existing reprodudion systems. One sudi unprovement that 

can be obtained In accordance with the invention is to use adaptive 
filters to compensate for practical deficiencies in loudspeaker 

loeatim relative to a listener. In the case of in-car 
entertainment systems for example, it is very difficult to locate the loudspeakers 
symmetrically with respect to the listener and thus produce the stereophomc 
illusion initially perceived by the recording engineer. The same may be true of 
many domestic listening environments. With appropriately designed filters, 
however, the loudspeakers can effectively be shifted to Virtual locations" which 
apparently seat the listener in die (^timal location. 

Th- o^ometrv used in computer simulations of this approach in 

I lorfre'St "'t'h'e f irst*5pect of the ^ V'^'^-jrilrand S2 to 

•> The matrix C(2) which relates the signals output from sources ?! 

miciophoJS^i and (pUced at the location of the listeners ears) is given by 













2-^22 


_ r2\ 





(16) 



where these transfer fundions are the digital versions of *e continu^^^^ 
transfer function relating the pressure at » P°«; ^.''"^ SvS by 
acceleration of a point monopole source. Thus the delays A,„| are given oy 

« round {f, r„i/c^ 

where r«, is the distance between Oie m th sourceand "jf 
the sampling frequency and c, is d« sound speed. ^^^'J^ 
spedfies the relationship between the recorded signals »«d Ae d«ired si^^ 
seSted by assuming that we have certain "virtual source »ocatK>ns ^^^^ 
« S^ted in FigiJe 2. This ensures that the desired signals di(n) and are 
thosethat would be produced by sources Viand V2. Thus 



A(z) = 



iii 



"11 «« 
■ "21 "22 



(18) 
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where z-Amarf is a modelling deUy and the delays Ymi given by 

= roundCf, u„d/ ^^^^ 

where u„, is the distance between the m'th virtual soura location axui ti^"*^ 
liophote. In aU the simulations that foUow, the sample rate /.was cho^ H, 
. ^ r>A^ ft - 1 0 ke/m^ for simpUdty. The number of 

Se^fi^^lt^S.: ™^hS) wlf i;^!;^ Chosen S be V and^e Biters 
^ SS^eS^adaptively using pseudorandom sequences x,(n) and x,M. 

in the arrangement shown in Figure ^Jtor^e^uxpo^^^ • 
sources SI a^d S2 are not only placed '^^fy'^^ ^ resulting 
but are also inverted relative to VI jnd i!^' HC^) has been designed 
impulse response of the sy'*«'" "^J w = I ^e corresponding 

using the algorithm of equation (14) with ^^'^'^^^-^J^^f i^ f dearly 
filter impulse responses are illt^^ated m ^^'^*„^^;,Jf^S^^^^^ i 
demonstrate that once signals mput to S^^nd S2 are premre y 
ON in the figures) the relative time of "J^^PfeT^ v\ and 

become equivalent to those * ^ rH^Sm^te Tht s'Sv^ fr^^^^ that 
V2. Note that for miaophone ^^^^P^f ''g^ before that from SI; 
from SI whilst for microphone 2 '•^^^ = OFF in 

a situation which reverses that f ^^^'^^^^^ ^J*S«t^a^i by the frequency 
the figure) The effectiveness of f ^J^^^^^'' ^'J^K how the relative 
response plots shown m Figure 5 ^^^^se m parat indicates the 

phfse beJween the signals arnvmg ^t'^^^J^^/J^S two microphones, 
difference in travel time between the ^2"*«*°St^t^^*Th^^ in 
Also note that the magnitude ^^'P^'^'^.'^'^^'^^S^^ Of course in this 
distance between the virtual sources the imoophor^ 
simple iUustration, no reverberation .s *^J"°**' """^ 

be acc=ounted for in the implementation of a real system. 

in accordance with a second ;^'>«'^^ .^^^^^ hst^^^^ wher. the 

such a technique could also find ^^P^'^^^^^^^^Zzd sources could 
loudspeakers for reproduction "^P'*«tJ^^ri)e^«n the real sources, 
be pl^«d in order to efff ^^^^ly inoease Ae s^a^^ ^^^^ 
Figures 7, 8 and 9 show the resuhs o «/»^2'^7~^^urces iUustrated in 
this approach with the geometry of *«/«f*.*"f' VfVr„uises applied via the 
Figure 6. Figure 7 shows "J'^^p^^S'Sf filt^matrix H(z). 
si^ xi(n) and are made different by ^'"^J;^^^ increased to be 
ln>rfcldar, the time betwe«i arnvaU at a^ i^^ 

consistent with the ^o^^^on of the vir tuad^^^ p.^, ^^,^^3 

filters necessary to accomphsh thw are H(z)C(r). These again 

virtual source locations. 

23 The production of mulHple stereophonic images 

In accordance with a third ^^P^f JSISS'mter design philosophy 

yet another possibiUty that emerges ^^^^^nnels of a conventional 

Outlined above, is that of 0P«"»"S °"JVde7Jirtual source locations for two 
stereophonic recording in order to produce ideal virm 
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independent listening positions. Such an approach does however require the 
use of at least four loudspeakers. As an illustration of this possibility, here we 
present the results of some computer simulations (first presented by Orduna- 
Bustamante et al [23]) where the source/microphone arrangement of Figure 10 is 
considered. Also, for the purposes of illustration, we will deal with the "cross- 
talk cancellation" case whae, for example, we wish to reproduce two recorded 
signals Zi(fi) and X2(n) at Ml and M2 and at M3 and MA. In this case, the matrix 
C(z) takes die form 



'(r^»/fn) (r^«/ri2) (r^»/r,3) (r^^^/ru) 

{X^ti/rn) (.r^^/f22> (r^»/r23) (r^24/r24) 

(r^31/r3i) (z-A32/r32) (r^33/r33) i^^/rii) 

Jr^UUx) (rA42/r42) (r^3/r43) (z^Vr44) . 

and the matrix A(z) which defines the desired signals is given by 



(20) 



A(z) = z-^mod 



1 0" 

0 1 

1 0 
LO 1. 



(21) 



where is a modelling delay. In the simulations presented here, the 4x2 
matrix H(z) was comprised of eight HR fillers each having a number of 
coefficients / = 128. The modelling delay ti„od was chosen to be 96 samples. 

The resulting impulse responses of the system consisting of H(z) C(z) are 
illustrated in Figure 11 and the corresponding frequency response functions are 
shown in Figure 12. A good degree of cross-talk cancellation is clearly evident, 
although the system becomes less effective at certain frequences where the 
inverse filters have difficulty in modelling what amount to poles of the inverse 
system (see Nelson et al {131 for a fuller discussion). The use of IIR fdters may 
provide a solution to this problem and a preliminary investigation is reported m 
reference 124). The use of HR filters in the single channel case has also been 
examined by Greenfield and Hawksford [251. This problem is more cleariy 
illustrated in the simple case of the reproduction of signals x-iin) and X2(n) at Ml 
and M2 respectively. The results of implementing such a system are illustrated 
in Figures 13 and 14. Effective cross-talk cancellation is produced (except at 
certain frequencies) at microphones 1 and 2, whereas the signals at microphones 
3 and 4 are clearly degraded by the operation of the system, as one would expect. 
Finally, the potential advantages in using more sources than microphones m the 
sense described in Section 2.3 is illustrated by the results presented in Figures 15 
and 16. Here, the 4 x 2 matrix H(z) has been designed to ensure reproduction of 
z,(«) and X2(n) at microphones 1 and 2. Considerably improved cross-talk 
cancellation can be seen to have been achieved. 
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SPARSELY UPDATED FILTERS FOR ADAPTIVE DIGITAL PROCESSING OF AUDIO 
SIGNALS 

A sparse update strategy will now be described that allows the 
implementation of adaptive filters at high sampling rates using 
existing DSP technology. This technique has the important property 
that the processing time per sampling period spent in filter update 
operations is independent of the filter length. Code will be given 
for both the LMS and the filtered-x LMS algorithms and description 
of their practical use for loudspeaker equalisation. 



The siiiglc^lutnnel system idaitification pn*lcm is shown m Figuxe 17Cai>.A vector 
5: ^ t«y«ifi g thft most Mcent samples of the inimt signal is defined as 

sunilady. a vector oomaining the coe€Eigf«ts of a (non-staiumaiy) non-feciirsive digital filter is 
defined as 

g- =[ go(n) giin) ^ go-tin) ]^ 
Tlie dot piDdua of the input vector with the coeffic^ 

This is leqaixed to approximate d« desired signal (the output ftom the system under 
idemification), in a vraty that inimimzcs the variance of the cn^ 

The LMS algorithm [31 pofbnmaswhasdc gndto 
famnla 

which can be shown m converge on die mean to the cxaa least me^ 
problem provided die adaptation caie is chosen to satisfy the fi^ 

\i<yLol 

The implementation of this aigocidun in a microprocessor requires 4L aritfmiedc operations 
(2L to perform the filtering and 2L to update the filter). M high samp to^ 
very taadng arid heavy restiicdons nam be imposed oh the number crfcoefBd 
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The system equaiiiaiioaimbleni is shown in Figute 17lbXInttaiscase.attewfi]tErwiih 

h.=[*o(n) ^ Aw(n) ]^ 
acts (m die input signal ID produce the sigi^ 

which after inuisniissiott through dte system produces an output that nuiumizBS the variance 
of a new error signal defined as 

The ooeffidents of die opdmal filter can be searohed using die filtend-x LMS algorithm [9 1 

hi.^-h.-i-)i£i.r«; 
whcro the vector r^, containing recent san^les of the reference signal 

is used instead of the vec»>r of input signals as in die LMS algocidim. This increases die 
ariduDBdc work by another 2L operadons (assuming diat bodi filters g and h have die same 
lengdi)> The arithnaedc woric thus increases to 6L ope n ti on s, 

SPARSE UPDATE IMPLEMENTATION OF THE LMS ALGORnHM. 

One way of cedudng the operation ccmnt re quir ed to update die filter ooeffidents 
using die LMS algoridun is to devise a criterion to select die abs(4ute imnimum number of 
operations that are stiU necessaiy to niaintain the convergence properties of the algotidmt The 
fTwnfftiimi work that can be done is to update only one filter coefficient per sampling period. 
This can be iiapiemeotBd by perfomiing die fonowing operadons at 

gA(«)«g*Ot)+Jie«x*^, (update current filter tap) 
k^lk-^iyaodL; Gncreoient tap counter) 

^i^iere ik is a counter set initially to h=0 that nms dicnlady along the vecUH' of filter 
coefficients. Note dial because n and t are always increriiented by I every 
k wraps around to zero), the input sample that is used for the update is exactly the same 
during die L cycles diat it takes to perform one pass along die whole filter. This observation 
\^iu\< to die fiiUowing aUeroatxve, but exacdy equivalent, versioo of the spmt update version 
of die LMS algoridun 

if I: s 0 then a s (ix. (stoic, and pre-multiply. input sample) 

gkM s gtOi) -i-ae« (update current filter op) 

its(t-i-l)nnod£ GncrementnpoouniBi) 

By following any of these procedures, it takes L proce s sing cydes to update the whole filter, 
but the operation otumt is reduced to 2L (basically those invoked in the actual filtering). 
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SPARSE UPDATE IMPLEMENTATlONpF THE FILTERED-X LMS ALGORITHM. 

» A sparse update itnpleoaentaiioii of the filtered-x LMS algorithm presents the 
additional challenge of having to calculate the filfcred rcfinence signal in a way which is 
compatible with the sparse update of the fdter coefficients. Interesdngly enough the 
calculation of the filtered reference signal can be pofonned also on a sin^e-tap bans as 
follows 



Note that the calculation of the next sample of the filtered reference signal starts L cydes in 
advance. To this end, the coefiBdents of the reference filter g are accessed in reverse, as 
shown, and die filtering tnakes use of the iDOSt rsoent inpm sample at ev^ 
were, gets old by itself as times goes by). Tbe operation count reduces again to just over 2L as 
in the qiarse update implementation of the LMS algortdun. 



LOUDSPEAKER EQUALIZATION USING SPARSE UPDATE ADAPTIVE FILTERS. 

Hgnre 18 shows die tequency response fonction of a loudspeaker in an anechoic 
chamber equaliaed to obtain a flat mHgmtnrie response and a linear phase response. The 
processing was peifoniMd in floatiiig poim aridunedc nang a Texas 

processor. The sampling frequency was set tD/^32 kHz and the filter length to L=A%. The 
impulse response fimctiao of dK^stem was first identified using die sparse update version of 
die LMS algDtidmi. The system was larerequaliaed using die sparse update im^ 
die filteied-x LMS algoridim. (A fuU update could only be pos^le at/slfikHz oc L=aA.) 



ifit^OdienrarO 

ifl: = L-l thcna = ^u• 
l:=(it+l)modL 



(dear filtered-x accumulator) 
(accumulate current product) 
(when done, store filiered-x) 
(increment tap counter) 



(update current filter U4>) 
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CLAIMS 

1. A sound reproduction system comprising a plurality of 
loudspeakers (SI, $2) for positioning at positions spaced from a 
listener location (Ml, M2), loudspeaker drive means (H) for driving 
the loudspeakers (SI, S2) in response to a plurality of channels of 
a sound recording (x), the loudspeaker drive means comprising a 
filter means (H), the sound recording (x) being suitable for playing 
normally through a plurality of speakers normally positioned 
generally symmetrically with respect to a listener location, 
characterised in that the filter characteristics of the filter means 
(H) are so chosen as to create, in use of the reproduction system, 
at the listener location (Ml, M2) a local sound field which is 
substantially equivalent to the local field that would result from 
playing the sound recording (x) through a plurality of loudspeakers 
driven without said filter means (H) and positioned at virtual 
source locations (VI, V2), at least one of the virtual source 
locations <V1, V2) being displaced substantially from the actual 
location (SI, S2) of the corresponding loudspeaker. 

2. A sound reproduction system as claimed in claim 1 characterised 
in that each of said virtual source locations (VI, V2) is displaced 
from the actual location (SI, S2) of the corresponding loudspeaker. 

3. A sound reproduction systan as claimed in claim 1 or claim 2 
characterised In that the actual locations (SI, S2) of the 
loudspeakers, during use of the sound reproduction system, are 
predetermined positions that are asymmetric with respect to the 
listener location (Ml, M2), whereas said virtual source locations 
(VI, V2) created by the filter means (H) are substantially 
symmetrically positioned with respect to the listener location (Ml, 
M2). 

4. A sound reproduction system as claimed in claim 2 characterised 
in that the actual locations (SI, S2) of the loudspeakers, during 
use of the reproduction system, are predetermined positions that ar 
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more narrowly spaced from each other than the spacing of the virtual 
sources (VI, V2) created by the filter means (H). 

5. A sound reproduction system comprising four loudspeakers (SI, 
S2, -Sa, S4) which are arranged in use at spaced-apart positions to 
create a sound field, first and second predetermined listener 
locations (Ml, M2), (M3, M4) within the sound field, characterised 
in that the loudspeakers are driven through filter means (H) by two 
channels (Xp Xg) of a sound recording, the filter characteristics 
of the filter means (H) being so chosen as to create at both the 
first and second listener locations (Ml, M2), (M3, M4) a respective 
local portion of the sound field which is substantially the same as 
the sound field portion that would be produced at that respective 
location (Ml, M2), (M3, M4) by playing the unfiltered channels of 
the sound recording (x^, Xg) through a pair of loudspeakers 
positioned symmetrically with respect to the respective location. 

6. An adaptive audio system comprising an adaptive filter having a 
plurality of alterable filter coefficients, a processor implementing 
an algorithm which cyclically performs in turn a filtering operation 
utilising the filter, and a filter updating operation to update the 
filter coefficients, characterised in that the algorithm is arranged 
in a cycle thereof to adjust (if necessary) only a limited number of 
the filter coefficients before performing a filtering operation. 

7. An adaptive audio system as claimed in claim 6 characterised In 
that only one coefficient of the filter is adjusted in each cycle. 

8. An adaptive audio system as claimed in claim 6 or claim 7 in 
which the algorithn is the LMS algorithm. 

9. An adaptive audio system as claimed in claim 6 or claim 7 
characterised in that the algorithm is the filtered-x LMS algorithm. 

10. An adaptive audio system as claimed in claim 9 characterised in 
that the filtered reference signal employed is calculated on a 
single-tap basis. 
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Figure 8 
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